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O ■ Abstract 

(N 

5h ■ We study spatio-temporal intermittency (STI) in a system of coupled sine circle maps. The phase 

"^ , diagram of the system shows parameter regimes with STI of both the directed percolation (DP) 

pg I and non-DP class. STI with synchronized laminar behaviour belongs to the DP class. The regimes 

of non-DP behaviour show spatial intermittency (SI), where the temporal behaviour of both the 

Q 

r J . laminar and burst regions is regular, and the distribution of laminar lengths scales as a power law. 

C I The regular temporal behaviour for the bursts seen in these regimes of spatial intermittency can 

be periodic or quasi-periodic, but the laminar length distributions scale with the same power-law, 
which is distinct from the DP case. STI with traveling wave (TW) laminar states also appears 

\^ , in the phase diagram. Soliton-like structures appear in this regime. These are responsible for 

vn 

^D I cross-overs with accompanying non-universal exponents. The soliton lifetime distributions show 

^^ ■ power law scaling in regimes of long average soliton life-times, but peak at characteristic scales 



o 
> 



with a power-law tail in regimes of short average soliton life-times. The signatures of each type 
of intermittent behaviour can be found in the dynamical characterisers of the system viz. the 
eigenvalues of the stability matrix. We discuss the implications of our results for behaviour seen 



S^ . in other systems which exhibit spatio-temporal intermittency. 
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I. INTRODUCTION 

The phenomena of spatiotemporal intermittency (STI), wherein laminar states which 
exhibit regular temporal behaviour co-exist in space and time with burst states of irregular 
dynamics, is ubiquitous in both natural and experimental systems. Such behaviour has been 
seen in experiments on convection ij, |2] , counterrotating Taylor-Couette flows |3|J , oscillating 
ferro-fluidic spikes [4] , experimental and numerical studies of rheological fluids J5|, |6[ , and 
in experiments on hydrodynamic columns |^. In theoretical studies, STI has been seen 
in PDEs such as the damped Kuramoto-Sivashinsky equation J8| and the one-dimensional 
Ginzburg Landau equation [9|, in coupled map lattices |lO| such as the Chate-Manneville 



ll|, the inhomogeneously coupled logistic map lattice |12], and in cellular automata 



CML 
studies. 

A variety of scaling laws have been observed in these systems. However, there are no 
definite conclusions about their universal behaviour. The type of spatiotemporal intermit- 
tency in which a laminar site becomes active (turbulent) only if at least one of its neighbour 
is active, has been conjectured to belong to the directed percolation (DP) universality class 
[13| . The dry state or the absorbing state in DP is identified with the laminar state in STI, 
and the wet state of DP corresponds to the active state in STI, with time as the directed 
axis. However, a CML specially designed to exhibit STI by Chate and Manneville showed 
critical exponents significantly different from the DP universality class [llj. This led to a 
long debate in the literature [iJ, ll^ lla, lllf . It was concluded that the presence of coherent 
structures, called solitons, were responsible for spoiling the analogy with DP. The nature 
of the transition to STI and the identification of the universality classes of STI is still an 
unresolved issue, and is a topic of current interest. 



Earlier studies of the diffusively coup 
which were completely free of solitons 



ed sine circle map lattice showed regimes of STI 
171 llq . Two types of STI were seen along the 
bifurcation boundaries of the bifurcation from the synchronized solution. The first type of 
STI showed an entire set of static and dynamic scaling exponents which matched with the 
DP exponents, and therefore was seen to belong convincingly to the DP class. The other 
type of intermittency, where both the laminar and burst regions showed regular temporal 
dynamics, was called spatial intermittency (SI). The laminar length distribution for this case 
showed characteristic power-law behaviour with its own characteristic exponent ( = 1.1. 



This kind of behaviour has been observed in the sine circle map lattice as well as in the 
inhomogenously coupled logistic map lattice. In the case of the sine circle map lattice, both 
types of intermittency, viz. STI of the DP class, and SI which does not belong to the DP 
class, were seen in different regions of the phase diagram. Moreover, distinct signatures of 
the two types of behaviour were picked up by the dynamical characterisers of the system, 
i.e. the eigenvalues of the stability matrix. The eigenvalue spectrum was continuous in the 
DP regime, but exhibited the presence of gaps in the SI regime. 

Different types of behaviour are seen within the SI class itself. The laminar state is 
synchronized in nature, but the burst state can be periodic or quasi-periodic in its dynamical 
behaviour. The periodic burst states can have different temporal periods. Burst states of the 
traveling wave type are observed at several points in the phase diagram. The distribution 
of laminar lengths shows power-law scaling in both cases with the same exponent. 

The SI regimes lie close to the bifurcation boundaries of the synchronized solutions. 
The SI with traveling wave (TW) bursts bifurcates further via tangent-period doubling 
bifurcations, to STI with TW laminar states and turbulent bursts. This type of STI is 
contaminated with coherent structures similar to the solitons that spoil the DP regime in 
the Chate-Manneville CML. The solitons induce cross-over behaviour and the exponents 
in this regime take non-universal values. The distribution of soliton lifetimes shows two 
characteristic regimes. In the short soliton lifetime regime, the distribution shows a peak 
which indicates the presence of a characteristic time scale, and has a power-law tail. In 
the longer soliton lifetime regime, the distribution has no characteristic scale and shows 
pure power-law behaviour. The solitons in this regime also change the order of the phase 
transition in the system. 

The dynamical characterisers of the system show signatures of the different types of 
temporal behaviour of the burst states. As mentioned earlier, the eigen-value distribution 
of the stability matrix is gapless for the STI with DP exponents, whereas distinct gaps are 
seen in the distribution for the SI class. The number of gaps in the eigenvalue distribution 
of SI as a function of bin size shows power-law behaviour. However, the scaling exponent 
is different for SI with quasi-periodic bursts and SI with periodic bursts. We discuss the 
implications of our results for behaviour seen in other systems which exhibit spatiotemporal 
intermittency. 

The organisation of this paper is as follows. Section II gives the details of the model and 
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FIG. 1: shows the phase diagram for the coupled sine circle map lattice evolved using random 
initial conditions. The spatiotemporally synchronized solutions are represented by dots. The points 
at which DP exponents have been obtained are marked by diamonds (O). At points marked with 
triangles (A), SI with quasiperiodic bursts is seen. SI with TW bursts is seen at points marked 
by crosses (x) and SI with period-5 bursts are seen at points marked with asterisks (*). STI with 
TW laminar states and solitons is seen at the points marked by boxes (□). The cluster solutions 
are marked with plus (+) signs . 

the phase diagram obtained. The two universality classes of spatiotemporal intermittency 
seen in this system, as well as the variations within the SI class are discussed in section III. 
Section IV explains the role played by the solitons in inducing a cross-over behaviour in STI 
with traveling wave laminar state. The signatures of each type of intermittent behaviour is 
seen in the dynamical characterisers of the system. This has been discussed in section V. 
We conclude with a discussion of these results and their implications for other systems. 

II. THE MODEL AND THE PHASE DIAGRAM 

The coupled sine circle map lattice is defined by the evolution equations 



X 



t+i 



^)f{x^ + ^[f{4-i) + f{xli)] (modi) 



(1) 



where i 



N and t are the discrete site and time indices respectively with A^ being 



the size of the system, and e being the strength of the coupling between the site i and its 



two nearest neighbours. The local on-site map, f{x) is the sine circle map defined as 

K 

f(x) =x + n- — sin(27rx) (2) 

ZTT 

Here, K is the strength of the nonlinearity and Q is the winding number of the single 
sine circle map in the absence of the nonlinearity. The coupled sine circle map lattice has 
been known to model the mode-locking behaviour 'j^ seen commonly in coupled oscillators, 
Josephson Junction arrays, etc, and is also found to be amenable to analytical studies J2i|. 
The phase diagram of this system is highly sensitive to initial conditions due to the presence 
of many degrees of freedom. Studies of this model for several classes of initial conditions 



m 



20|. 



have yielded rich phase diagrams with many distinct types of attractors 

We study the system with random initial conditions. The system is updated syn- 
chronously with periodic boundary conditions in the parameter regime < Q < ^ and 
K = 1 (where the single circle map has temporal period 1 solutions in this regime); the 
coupling strength, e is varied from to 1. 

The phase diagram obtained using random initial conditions is shown in Figure ^ J21|. 
Spatially synchronized, temporally frozen solutions, where the variables Xi{t) take the value 
Xi{t) = 3;* = 2~: sin~^ ^^ for all i = 1, . . . ,N, and for all t, are marked by dots in Figure 
E These solutions are seen over a large section of the phase diagram and are stable against 
perturbations. Cluster solutions, in which Xi{t) = Xj{t) for all i,j belonging to a particular 
cluster, are identified by plus signs (+) in the phase diagram. Regimes of spatiotemporal 
intermittency of various kinds are seen near the bifurcation boundary of the synchronized 
solutions. The various types of STI seen are: 

i. STI of the type in which the laminar state is the synchronized fixed point x* defined 
earlier, and the turbulent state takes all other values other than x* in the [0, 1] interval, 
is seen at points marked with diamonds (O) in Figure H The space-time plot is shown 
in Figure |2l^a). This type of STI belongs to the directed percolation universality class. 

ii. STI with TW laminar state interspersed with turbulent bursts is seen at points marked 
with boxes (□) in Figure ^ The space-time plot of this type of solutions is shown 
in Figure Efb). Coherent structures traveling in space and time are seen in these 
solutions. Such structures have also been seen in the Chate - Manneville CML and 
have been called solitons in the literature. 
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FIG. 2: shows the space time plots of the different types of STI seen in the phase diagram. 
The lattice index i is along the x-axis and the time index t is along the y-axis. The space time 
plots show (a) STI with synchronized laminar state interspersed with turbulent bursts seen at 
= 0.06, e = 0.7928. (b) STI with TW laminar state and turbulent bursts with solitons seen 
at = 0.037, e = 0.937. (c) SI with synchronized laminar state and quasi-periodic bursts seen 
at fi = 0.031, e = 0.42. (d) SI with synchronized laminar state and TW bursts observed at 
O = 0.019, e = 0.9616. 



iii. Spatial intermittency with synchronized laminar state and quasi-periodic bursts are 
seen at parameters marked with triangles (A) in the phase diagram. The space-time 
plot of this type of solution is shown in Figure |21^c). 

iv. Spatial intermittency with synchronized laminar state and traveling wave (TW) bursts 
are seen at points marked with crosses ( x ) in the phase diagram. The space-time plot 
is shown in Figure Efd). SI with synchronized laminar state and period-5 bursts are 
seen at points marked with asterisks (*) in the phase diagram. 

The identification of the universality classes of the different types of intermittency seen 
in this system has been partially carried out earlier. STI with synchronized laminar states 
and turbulent bursts has been clearly established to belong to the directed percolation (DP) 



class 







18|. However, the other types of intermittency seen in the parameter space do not 



belong to the DP class. We analyse these in the next section. 

III. THE UNIVERSALITY CLASSES IN THE SYSTEM 

It is interesting to note that the system under study exhibits spatiotemporal intermittency 
belonging to distinct universality class at different values of the parameters. The two distinct 
classes obtained so far are regimes of STI which belong to the DP class and regimes of SI 
which do not belong to the DP class. We discuss each of these in further detail. 

A. STI of the DP type 

It has been shown convincingly that STI with synchronized laminar state interspersed 
with turbulent bursts (seen at points marked with O in Figure lU) belongs to the DP univer- 

nn 

sality class [17|, |18|| . The infective dynamics of the turbulent bursts wherein the turbulent 
site either infects the adjacent laminar sites, or dies down to the laminar state, is similar to 
the behaviour seen in directed percolation J22| . Since no spontaneous creation of turbulent 
bursts takes place, the laminar state forms the absorbing state and time acts as the directed 
axis. The entire set of static and dynamic scaling exponents obtained in this parameter 
regime match with the DP exponents. The exponents obtained, after averaging over 10^ 
initial conditions, at three such parameter values are listed in Table HI The complete set of 



exponents and their definitions have been reported in |l7l . |l8[ . The distribution of laminar 
lengths also shows a scaling behaviour of the form, P{1) ~ Z^*', with an associated exponent, 
( ~ 1.67. The laminar length distribution obtained, after averaging over 50 initial condi- 
tions, at fi = 0.06, e = 0.7928 has been plotted in FigurelSfa). The size of the lattice studied 
was 10"^. 

A clean set of DP exponents is obtained for the STI with synchronized laminar state 
seen in this parameter regime as these regimes are completely free from the presence of 
coherent soliton-like structures which could bring in long-range correlations in the system 
and thereby spoil the DP behaviour. In fact, in the case of the STI with synchronized 
laminar state, no solitons have been observed for this model in the range of parameters 
studied. However, the STI with traveling wave laminar states seen at the parameter values 
marked with boxes in the phase diagram, does show the presence of solitons as seen in the 
space-time plot of Fig|2fb). These solitons are responsible for non-universal exponents and 
cross-over behaviour in this regime. This behaviour is discussed in detail in section IV. In 
the remainder of the present section we will discuss the second universality class seen for 
the present model, viz. that of spatial intermittency. 



Static and dynamic scaling exponents for the STI of the DP class 
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TABLE I: The static and dynamic exponents obtained in the DP regime at three of the points 
marked with diamonds (O) in the phase diagram are shown in the above table. The universal 
DP exponents are listed in the last row. The exponents, /? and v have been calculated using the 
hyperscaling relations. See [l3| for definitions of the DP exponents. 






(b) 


a 


'""-■,. 




"■"T--^ 


-1 


'"""'""*"'*^ 


-2 


^"^^Sv^ 


r^ 


^^^■-.^ 


K-3 


N. " ,. 


3 


\ 


-i 


\ 


-i 


\ 


-y 


1 II i 



Qi 



Hi 



>10 



25 




FIG. 3: shows the log — log (base 10) plot of the laminar length distribution for (a) STI of the 
DP class obtained sA Q. = 0.06, e = 0.7928. The exponent C, is 1.681. (b) SI with quasi-periodic 
bursts obtained at fi = 0.04, e = 0.402. The exponent C, is 1.10. (c) SI with TW bursts obtained 
at r^ = 0.019, e = 0.9616. The exponent obtained is 1.08. 

B. Spatial Intermittency 



We now discuss the other type of intermittency, seen in this system viz. spatial intermit- 
tency. Spatial Intermittency is a distinct class of STI in which the temporal behaviour of 
both the laminar and burst states is regular. The infective dynamics characteristic of the 
DP class is absent here, and the burst states do not infect the laminar states even when they 
are nearest neighbours. Spatial intermittency is a long lived phenomenon and the spatially 
intermittent state persists for time scales which are much longer than the time scales on 



Spatial intermittency with quasi-periodic bursts 
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TABLE II: The table shows the laminar length distribution exponent, ^, calculated for SI with 
quasi-periodic bursts (marked by triangles (A) in Figure^. The exponent k is the exponent 
associated with the number of gaps, Ng{l) ~ l^'^ in the eigenvalue distribution, where / is the bin 
size chosen. The frequencies, lo inherent in the time series of the burst state are also listed and are 
of the form toi, 0^2, and toi + 102- 

which the STI states die down to a uniform laminar background. Two different types of 
spatial intermittency have been seen in this system. In both types of SI, the laminar state 
is the synchronized fixed point, x* defined earlier. However, the burst states are different 
and may be either quasi-periodic (marked by triangles in the phase diagram) or periodic in 
their temporal behaviour. 

1. SI with quasi-periodic bursts 

Spatial intermittency, in which the temporal behaviour of the burst states is quasi- 
periodic in nature, has been seen at points marked with triangles (A) in Figure^ The 
space-time plot has been shown in Figure |2fc). These burst states are non- infective in 
nature, i.e. the probability of the burst state infecting the nearby laminar state is zero. 
Therefore, the laminar states remain laminar forever. Hence, after an initial transient, 
the order parameter of the system, which is defined as the fraction of non-laminar sites 
in the lattice, is a constant. The time series of the burst states at different parameter 
values, at a typical burst site, was studied using power spectrum analysis. The power 
spectrum obtained at i7 = 0.058, e = 0.291 and Q = 0.0495, e = 0.3178 has been shown 
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FIG. 4: shows the power spectrum , |F(lij)P of the time series of the burst state seen at (a) 
n = 0.058, e = 0.291 and (b) n = 0.0495, e = 0.3178. The time series shown in (a) exhibits 
quasi-periodic behaviour and the time series shown in (b) is periodic in nature. 

in Figure m As can be seen from Figure Ufa), the peaks are seen at ui, UJ2, uji + U2 
and at mui + nuj2. This kind of behaviour is typical of a quasi-periodic state. Hence, 
we confirm that a quasi-periodic burst state is seen at Q = 0.058, e = 0.291. 

The laminar length distribution of this type of SI shows a scaling behaviour of the form, 
P{1) ~ l^'' with an associated exponent, ( ~ 1.1. The laminar length distribution at 
parameters {Q = 0.04, e = 0.402) has been plotted in Figure Efb). The values of ( 
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FIG. 5: shows the order parameter, m vs (e — ec) plotted for SI with synchronized laminar state 
and quasiperiodic bursts at J7 = 0.031 (A) and STI of the DP class at il = 0.06 (O) in Fig. (a). 
Fig (b) shows m vs (e — Cc) plotted for SI with TW bursts at Jl = 0.042. SI shows a first order 
transition whereas DP class shows a second order transition. The data collected is for a 1000 site 
lattice and is averaged over 1000 initial conditions. 
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Spatial intermittency with 


periodic bursts 




n 


e 


c 


K 


CJ 


0.019 


0.9616 


1.08 ± 0.04 


1.18±0.04 


0.5 


0.025 


0.9496 


1.08 ± 0.02 


1.12±0.02 


0.5 


0.037 


0.9254 


1.17 ± 0.02 


1.07±0.01 


0.5 


0.042 


0.9148 


1.13 ± 0.02 


1.10±0.04 


0.5 


0.047 


0.3360 


1.13 ± 0.02 


1.02±0.01 


0.2, 0.4 


0.0495 


0.3178 


1.15 ± 0.04 


1.03±0.02 


0.2, 0.4 


0.054 


0.2936 


1.17 ± 0.03 


1.02±0.02 


0.2, 0.4 



TABLE III: The table shows the laminar length distribution exponent, C calculated for SI with 
periodic bursts (marked by crosses (x) and asterisks (*) in Figure ^). The exponent k is the 
exponent associated with the number of gaps, Ng{l) ~ l^'^ in the eigenvalue distribution, where 
/ is the bin size chosen. The frequencies, uj inherent in the time series of the burst state are also 
listed. 

obtained for this type of SI at different values of {Q, e) have been listed in Table |nj 
The scaling exponent, ( obtained for this type of SI is clearly different from that of 
the DP class {(dp = 1-67). 

Secondly, the transition to SI from a completely synchronized state is a first order 
transition unlike the transition to STI belonging to the directed percolation class which 
shows a second order transition. This can be seen in Figure E^a) in which the order 
parameter of the system m, which is defined as the fraction of burst states in the lattice, 
has been plotted as a function of the coupling strength, e. The order parameter, m 
increases continuously with e in the case of STI of the DP class, signalling a second 
order transition, whereas m shows a sharp jump with e in the case of SI with quasi- 
periodic bursts, indicating that a first order transition takes place in the case of SI. 

2. SI with periodic bursts 

Spatial intermittency with periodic bursts forms the second class of SI. Two distinct 
burst periods have been observed in the phase diagram. Bursts of period 5 have been 
seen at the points marked by asterisks (*) in the phase diagram. 

12 



Figure Efb) shows that peaks are seen in the power spectrum of the burst state time 
series at u; = 0.2 and higher harmonics. This confirms that the burst states have 
period 5 at fi = 0.0495, e = 0.3178, which is one of the points marked by asterisks in 
the phase diagram. 

Bursts of spatial period two, temporal period two, of the traveling wave (TW) type are 
seen at the points marked by crosses (x) in the phase diagram. The laminar state in 
both cases is the spatiotemporally synchronized fixed point, x*. The space-time plot 
of these SI with TW burst solutions at 1] = 0.019, e = 0.9616 is shown in Figure Efd). 
The bursts are non- infective in nature in this type of SI as well. The scaling exponent, 
( associated with the laminar length distribution at different parameter values in this 
regime have been listed in Table IIIII for both TW and period 5 bursts. The laminar 
length distribution exponent obtained in this regime is ( ~ 1.1. The transition to SI 
with TW burst state from a spatiotemporally synchronized state is also a first order 
transition as has been shown by the abrupt jump in the order parameter, m with 
change in the coupling strength, e (Figure E^b)). 

It is thus clear that SI does not belong to the DP universality class. The scaling 
exponent ( = 1.1 for laminar lengths for the SI is distinctly different from the DP 
exponent ( = 1.67. We note, however, that the nature of the bursts, viz. periodic or 
quasi-periodic, has no effect on the value of the exponent (. We hence conclude that SI 
with periodic as well as quasi-periodic bursts belong to the same class. A similar value 
of the laminar length distribution exponent {( ~ 1.1) has been reported for spatial 
intermittency in the inhomogenuously coupled logistic map lattice 12|- Thus spatial 



intermittency appears to constitute a distinct universality class of the non-DP type. 

Therefore, two distinct universality classes of spatiotemporal intermittency, viz. directed 
percolation and spatial intermittency, are seen in the coupled sine circle map lattice in 
different regions of the parameter space. The reasons for the appearance of these two 
distinct classes may lie in the long-range correlations in the system at different parameter 
values. 
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Bifurcations from SI with TW bursts 



ft 



Eigenvalues Type of bifurcation 



0.0100 0.982 1.685, -1.618 

0.0210 0.960 1.361, -1.222 

0.0305 0.943 1.752, -1.535 

0.0410 0.920 1.623, -1.309 



TP 
TP 
TP 
TP 



TABLE IV: shows ttie largest positive and largest negative eigenvalues observed when SI with 
TW bursts (marked with crosses in the phase diagram) bifurcates to STI with TW laminar states 
(marked with boxes). The solution changes through a tangent-period doubling (TP) bifurcation. 

IV. THE ROLE OF SOLITONS IN STI WITH TW LAMINAR STATE 

As mentioned in Section III. A, in addition to spatiotemporal intermittency with synchro- 
nized laminar states, the phase diagram of our model also shows spatiotemporal intermit- 
tency with TW laminar states and turbulent bursts at the points marked by boxes in the 
phase diagram. The lattice dies down to the absorbing TW laminar state from random ini- 
tial conditions asymptotically. The STI with TW laminar states seen in this model appears 
as a result of a tangent-period doubling bifurcation from the SI with TW bursts as can be 
seen from Table IIVI 

Immediately after the bifurcation, apart from turbulent states, coherent structures, which 



Laminar length exponents in the solitonic regime 





Q = 


= 0.035 




n = 


= 0.037 


e 




c 


e 




c 


0.933 




1.53 ± 0.01 


0.930 




1.50 ± 0.02 


0.943 




1.40 ± 0.01 


0.937 




1.31 ± 0.05 


0.950 




1.17 ± 0.01 


0.950 




1.17 ±0.01 


0.962 




1.02 ± 0.01 


0.962 




1.02 ± 0.00 



TABLE V: shows the laminar length distribution exponent, ^ obtained for different values of the 
coupling strength, e and for different il's in the STI with TW laminar state and turbulent bursts 
regime. 
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100000 




FIG. 6: shows the log-log (base 10) plot of escape time, r plotted as a function of the lattice size, 
L at n = 0.037 and at e = 0.937 (O), 0.938 (□), and at e = 0.939 (A). The data is averaged over 
1000 initial conditions. 

have been called solitons, are seen in the TW laminar background. These structures have 
been marked in the space-time plot of this type of STI in Figure |2fb). The solitons travel 
through the lattice with a velocity, v = 1/t' such that for a right moving soliton, x* = 
xl'l^l , and X* = x*lj for a left moving soliton. Here, i and t are the site and time indices 
respectively. In this model, the left and right moving solitons occur in pairs and hence, 
they annihilate each other. When these solitons collide, they either die down to the TW 
laminar state or give rise to turbulent bursts. Such coherent structures have been seen 
earlier in the Chate Manneville CML jl4| , where these solitons were responsible for spoiling 
the DP behaviour and were even capable of changing the order of the transition. In the 
case of the sine circle map CML as well, the solitons seen in the STI with TW laminar 
state are responsible for non-universal exponents with values varying from 1.0 to 1.5 in the 
distribution of laminar lengths (Table IVJl . Additionally, the escape time, r which is defined 
as the time taken for the lattice to relax to a completely laminar state, starting from random 
initial conditions, does not show power-law scaling as a function of L for any value of the 
parameters (See Fig. ^. 

The soliton life-times and velocities depend on the coupling strength e and fl. As the 
coupling strength, e increases, the velocity, v of the solitons is found to increase. Therefore, 
solitons with larger velocities collide earlier with each other, and hence have shorter lifetimes. 
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FIG. 7: (a) shows the log-log (baselO) plot of distribution of soliton lifetimes at 17 = 0.035 and 
at e = 0.933 (C = 1.53), e = 0.943 (C = 1.40), e = 0.95 (C = 1.14) and at e = 0.962 (C = 1.02). 
The exponent, ( depends on the soliton lifetimes, (b) shows the log-log (baselO) plot of the 
soliton lifetime distribution at fi = 0.035, e = 0.943 (diamonds); 0, = 0.037, e = 0.937 (pluses) 
and at r2 = 0.04, e = 0.9323 (boxes). The laminar length distribution exponent, ^ is -^ 1.3 at all 
these points. Hence, the soliton lifetimes collapse on to each other, (c) shows the soliton lifetime 
distribution seen in the short soliton lifetime regime at fi = 0.037, e = 0.962. A peak is seen at the 
characteristic timescale, Tc = 20. The distributions have been obtained after evolving a lattice of 
size 500 over 20, 000 time steps and 50 initial conditions. 

The distribution of soliton life-times shows a peak in the short life-time regime indicating 
the presence of a characteristic soliton life-time r^. However, the tail of the distribution falls 
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off with power-law behaviour with an exponent 2.84. For low values of e, where the soliton 
velocities are smaller, there is no peak or characteristic life-time in the distribution, and the 
entire distribution of soliton life-times scales as a power-law with exponent ~ 1.1 (See Fig. 
|7)) . It is seen that the exponent ( for the laminar lengths decreases as the life-times decrease 
and the turbulent spreading in the lattice decreases. 

The spreading dynamics in this type of STI was studied by introducing a cluster of 
turbulent seeds in a completely absorbing background. The two dynamic quantities (i) N{t) 
the fraction of turbulent sites in the lattice at a time t, and (ii) the survival probability, P{t), 
which is defined as the fraction of initial conditions at time t which show a non-zero number 
of active sites, were studied at Q = 0.035 and e = 0.933, 0.943, 0.95, and 0.962. These have 
been plotted in Figure IHa) and (b) respectively. It can be seen from the figure that the 
fraction of turbulent sites, N(t) at a given time, t decreases as the coupling strength, e is 
increased. We see a similar decrease in the fraction of initial conditions which survive, P{t) 
with increase in e. The data is averaged over 1000 initial conditions. 

Therefore, we can conclude that the extent of spreading in the lattice decreases as the 
lifetime of the soliton decreases (with increase in e). Since the distribution of laminar lengths 
is an indirect measure of the spreading in the lattice, we see that the varying average soliton 
lifetimes influence the distribution of laminar lengths. Therefore, the solitons seen in this 



Solit( 


on lifetimes in STI with TW laminar state 




Regime 


n 


e 


c 


Tc 


-^ max 


/" 


Long soliton 


0.035 


0.933 
0.943 


1.53 
1.40 


- 


19010 

2481 


1.14 
1.35 


lifetimes 


0.037 


0.930 
0.937 


1.50 
1.31 


- 


16961 

5782 


1.18 
1.31 


Short soliton 
lifetimes 


0.035 
0.037 


0.962 
0.962 


1.02 
1.02 


15 
20 


305 
413 


2.84 
2.88 



TABLE VI: shows the exponents obtained in the short and long soliton lifetime regimes. Here, C is 
the laminar length distribution exponent. A characteristic timescale, Tc, is seen in the short soliton 
lifetime regime. Tfnax is the largest soliton lifetime observed and /x is the exponent associated with 
the soliton lifetime distribution. 
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FIG. 8: shows the log-log plot of (a) the fraction of turbulent sites, N{t) plotted as a function of t 
at $7 = 0.035, and (b) the survival probability, P{t) plotted as a function of t at 17 = 0.035. The 
extent of spreading decreases with increasing coupling strength, e or equivalently with decreasing 
soliton lifetimes. 

regime are responsible for non-universal exponents here. Conversely, the soliton lifetime 
distributions have been plotted in Figure [3 (b) for parameters where the exponents ( for the 
laminar length distributions take similar values. The soliton lifetime distributions collapse 
over each other as expected. 

We note again that the STI with TW laminar state shows no soliton free regime, and 
the DP regime where the laminar state is the synchronized state is completely soliton free. 
Hence, no direct comparison of the exponents of the STI with synchronized laminar state 
and STI with TW laminar state is possible at present. 

V. DYNAMIC CHARACTERISERS 

It has been seen earlier that the signature of the DP and non-DP behaviour in this sys- 
tem can be seen in the dynamical characterisers of the system, specifically in the eigenvalues 
of the one-step stability matrix. Regimes of STI with DP behaviour exhibited a contin- 
uous eigenvalue spectrum, whereas regimes of spatial intermittency showed an eigenvalue 
spectrum with level repulsion, where distinct gaps were seen in the spectrum. 

The linear stability matrix of the evolution equation^at one time-step about the solution 
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FIG. 9: shows the eigenvalue distribution (bin size=0.005) for (a) STI belonging to the DP class at 
0, = 0.06, e = 0.7928, (b) Spatial intermittency with quasi-periodic bursts at Q = 0.04, e = 0.402, 
and (c) SI with TW bursts at Q = 0.025, e = 0.9496. A section of the eigenvalue distribution is 
magnified in the inset figures. Gaps are seen in the spatial intermittency eigenvalue distributions 
whereas the eigenvalue distribution for STI does not show any such gaps. 



of interest is given by the N x N dimensional matrix, M^ , given below 

^ esf'ixi) ej'{x\) ... ej'{x%) \ 

ej'{x\) esf'ixi) ej'ixl) ... 

e„f(4) ej'(x*) ... 



M 



N 



[enf'ix\) 







e„/'(a;^_i) ej'{x%) J 
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where, eg = 1 — e, e„ = e/2, and f'{xj) = 1 — K 008(271x1) . x\ is the state variable at site i 
at time t, and a lattice of A^ sites is considered. The diagonalisation of the stability matrix 
gives the N eigenvalues at time t. 

A. The eigenvalue distribution 



The eigenvalue distribution, P(A) of the matrix above, in all the cases studied here, have 
been obtained after averaging over 50 initial conditions for 1000 lattice sites. Figure IHl shows 
the eigenvalue distributions for STI belonging to DP class at the typical parameter value 
VL = 0.06, e = 0.7928 (a), SI with quasi-periodic bursts at the typical value Q = 0.04, e = 
0.402 (b), and SI with TW bursts at fi = 0.025, e = 0.9496 (c). The bin size chosen is 0.005. 
It is clear from the figure that the eigenvalue distribution of the DP class at this value of 
bin size is continuous whereas distinct gaps can be seen in the distribution for the spatial 
intermittency class for both quasi-periodic and periodic bursts. 

For the case of the SI, the number of vacant bins in the eigenvalue distribution, Ng(l) 
scales as a power-law Ng{l) ~ Z^" where / is the bin-size f Figure ITUj). However, the exponent, 
K depends on the inherent dynamics of the burst states. The exponent, n for SI with quasi- 
periodic bursts have been listed in Table |n] and the exponents for SI with periodic bursts 
are listed in Table IIIII . 

Within the SI class, the value of k is seen to be stable within each period for the periodic 
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FIG. 10: shows the number of vacant bins, Ng{l) plotted against the bin size, I on a log- log (base 
10) scale for (a) SI with quasi-periodic bursts, and (b) SI with periodic bursts. The exponent, k 
associated with Ng{l) ~ l^'^ is given in the figures. 
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FIG. 11: shows the power spectrum of the time series of the largest eigenvalue, Xm{t) at (a) 
il = 0.058, e = 0.291, where SI with quasi-periodic behaviour is seen, (b) Q = 0.06, e = 0.7928, 
where STI of the DP class is seen, and (c) il = 0.026, e = 0.948, where SI with TW bursts is seen. 
Three main frequencies wi, uj2 and ui +uj2 are seen in (a) whereas, a broadband spectrum is seen 
in the case of STI of the DP class as seen in (b). The spectrum of SI with TW shows a peak at 
0.5 as seen in (c). 

bursts (See Fig. ITO|) . In the quasi-periodic case, the natural frequencies of the dynamics are 
different at different values of the parameter and hence k values are different for different 
values of the parameters. It is also useful to track the temporal evolution of the largest 
eigenvalue to identify the signatures of the differences between these three cases. 

B. The temporal evolution of the largest eigenvalue 

The temporal evolution of the largest eigenvalue of the stability matrix, A^, with t contains 
information about the dynamical behaviour of the burst states. The time series of the largest 
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eigenvalue was obtained for the three cases: STI of the DP class, SI with quasi-periodic 
bursts, and SI with TW bursts. After the initial transient, Xm settles down to the natural 
periods of the burst states. The power spectrum, |F(ci;)p picks out the inherent frequencies 
in the system. This is evident from Figure ITD in which the power spectrum of the time series 
of Xm(t) has been plotted as a function of the frequencies. 

In the case of SI with QP bursts (Figure HT]^ a)), peaks are seen at ui, uj2, uji + uj2 and 
at muJi + nuj2 {m, n > 0) as is typical of a quasi-periodic behaviour. However, in the case 
of STI of the DP class , a broad band spectrum is obtained which implies that the burst 
states contains many frequencies. In the case of SI with periodic bursts, peaks are seen in 
the power spectrum aX uo = 0.5 for SI with TW bursts f Figure ITlT c)). indicating period-2 
temporal behaviour, whereas the peaks are seen at u; = 0.2, and 0.4 for the type of SI in 
which the temporal behaviour of the bursts is period-5. 

Thus we note that DP behaviour is associated with a broad-band spectrum for the power 
spectrum of the temporal evolution of the largest eigenvalue, as well as a gapless distribution 
of eigenvalues, whereas the SI or non-DP behaviour is associated with the characteristic 
power spectrum of the temporal nature of the burst states i.e. periodic or quasi-periodic 
behaviour, and distinct gaps in the eigenvalue distribution. 

VI. CONCLUSIONS 

Thus, spatiotemporal intermittency of several distinct types can be seen in different 
regions of the phase diagram of the coupled sine circle map lattice. STI is seen all along the 
bifurcation boundaries of bifurcations from the synchronized solutions. These bifurcations 
are of the tangent-tangent (TT) and tangent-period doubling (TP) type. The universal 
behaviour of the system as typified by the laminar length exponents is of two types- the 
directed percolation (DP) class and the non-DP class. STI with synchronized laminar states 
belongs convincingly to the DP class and can be seen after both TT and TP bifurcations 
from the synchronized state. This class of STI is remarkably free of the solitons which 
spoil the DP behaviour in other models such as the Chate Manneville CML. Other regimes 
of the phase diagram show spatial intermittency (SI) behaviour where the laminar regions 
show power-law scaling and are periodic in behaviour, and the burst states show temporally 
regular behaviour of the periodic and quasi-periodic type. This type of intermittency is 
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clearly not of the DP type and has earlier been seen in the inhomogenously coupled logistic 
map lattice. 

In addition to the two regimes above, we also see STI with traveling wave (TW) laminar 
state in some regions of the parameter space. This kind of STI arises as the result of a TP 
bifurcation from SI with synchronized laminar state and TW bursts. This type of STI is 
contaminated with solitons and hence shows non-universal exponents. The soliton lifetimes 
depend on the parameter values and their distributions show two characteristic regimes. In 
the first regime, where typical lifetimes are short, the distribution peaks at short lifetimes 
showing the presence of a characteristic soliton lifetime scale but has a power-law tail. 
In the second regime, where soliton lifetimes are typically larger, the distribution has no 
characteristic scale and shows power-law behaviour with an exponent in the range 1.1 — 1.2. 

The dynamic characterisers of the system, namely, the eigenvalues of the stability matrix, 
shows signatures of these distinct types of behaviour. The DP regime is characterised by 
a gapless eigenvalue distribution and a broadband power spectrum of the time series of 
the largest eigenvalue. For the SI case, i.e. the non-DP regime, distinct gaps are seen 
in the eigenvalue distribution, and the power spectrum of the temporal evolution of the 
largest eigenvalue is characteristic of periodic or quasi-periodic behaviour depending on the 
temporal nature of the burst states. 

The origin of the different types of universal behaviour in different parameter regimes 
appears to lie in the long range correlations in the system. These correlations in the system 
appear to change character in different regimes of parameter space leading to dynamic 
behaviour with associated exponents of the DP and non-DP types. In order to gain insight 
into the nature of the correlations in this system, and the way in which they change character 
in different parameter regimes, we plan to set up probabilistic cellular automata which 
exhibit similar regimes and to examine their associated spin Hamiltonians l23||. Absorbing 
phase transitions are seen in other CML-s j24l . 12 5L I26| and in pair contact processes [22, 



l28l |. Models of non-equilibrium wetting set up using contact processes with long-range 
interactions also show DP or non-DP behaviour depending on the activation rate at sites at 



the edges of inactive islands 



Similar ideas may apply to the behaviour in our model as 



well. We hope to explore this direction in future work. 
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